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ABSTRACT 

A through-flow (hub-to-shroud) truly inverse method is proposed and 

described in this paper. It uses as a design initial specification an imposition of mean 
swirl, i.e., radius times mean tangential velocity, given throughout the meridional 
section of the turbomachine. In the present implementation, it is assumed that the 
fluid is inviscid, incompressible and irrotational at inlet and the blades are supposed to 
have zero thickness. Only blade rows that impart to the fluid a constant work along the 
span will be considered. 

An application of this procedure to design the rotor of a mixed-flow pump will 

be described in detail. The strategy used to find a suitable mean swirl distribution and 

the other design inputs is also described. The final blade shape and pressure 
distributions on the blade surface are presented, showing that it is possible to obtain 

feasible designs using this technique. Another advantage of this technique is the fact 

that it does not require large amounts of CPU time. 


1 -INTRODUCTION AND LITERATURE SURVEY 

The large majority of pumps is designed by using very simple and rudimentary 
one-dimensional considerations concerning the velocity triangles, considerations 
which allow the calculation of the evolution of the blade angle along the passage, see 
for example [1] and [2], These methods are so easy that they can be carried out using 
only hand calculations and simple graphical processes. 

Although some pumps are still being calculated using hand calculations, the 

above one-dimensional procedure can and has been programmed as computer codes 
which are being used by the most important pump manufacturers. The results 
obtained with these techniques are the better the more radial the blade passage is. 
However, for mixed-flow pumps, the velocity triangles vary appreciably along the 

span, so that the above methodology is not good enough. In order to take into account 
this effect in some way, designers usually split the flow passage in several parts, 
applying the above considerations to each one. 

When the pump designs involve some responsibility, the above step is followed 
by a verification using a direct code which is run with the geometry arrived at 
previously. Before reaching the final design, several iterations following the above 
steps are usually required. 

This complete process can be time consuming, so that one is left wondering 

whether it could be improved. One possible way to achieve this could be by using 
inverse methods, enabling the achievement of the blade row layout in a more direct 
form, in one single step. Among inverse methods, two-dimensional techniques are the 
most frequently used and the ones that require less CPU time. These two-dimensional 
inverse methods can be classified into two main groups, according to the 
approximations used when looking for the solution to the blade design. The first sort of 
approximations gives rise to blade-to-blade methods since the calculations are done in 
the blade-to-blade plane. This kind of methods is popular among designers of axial 
turbomachinery, but it did not attract much attention among designers of radial 
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turbomachinery. The reason may be connected to the fact that the flow passages are 
more complicated and there are significant changes in radius. 

The other way of tackling the solution consists in working in the hub-to- 

shroud plane (i.e., in the meridional plane), giving rise to what we could call hub-to- 

shroud (or through-flow) methods. In this sort of approach the flow is supposed 
axisymmetric, an approximation that can be interpreted as giving the mean flow 

through the turbomachine. This sort of approach seems more popular among radial 

turbomachinery designers. In fact, one can find examples of this kind of procedure as 

far back as 1955, when reference [3] described some work applied to the redesign of a 
centrifugal compressor, using a hub-to-shroud method. In this instance the blade 
shape was kept fixed and given as input, while the shroud contour was altered and 

evolved as a result of the calculations. The centrifugal compressors described in [3] 
were built and tested, the experimental results being presented in [4]. The 
experimental results show that this technique produced significant improvements in 

the overall efficiency and peak pressure ratio. 

A different strategy was followed in reference [5] which, again, describes the 
application of an inverse technique to the design of centrifugal compressors. 

Contrary to the previous example the meridional geometry (hub and shroud contours) 
of the machine is supposed known and given as input. The other inputs consist of a 

suitable normal blade thickness and the desired velocity loading (difference in 
velocity across the blade) at hub, mid-span and shroud as a function of distance along 

the camberline. As a result of the calculations the blade shape was obtained. This 
paper presents some experimental evidence suggesting that the procedure gives 
reliable results when the flow is attached. 

The method proposed in the present work has got some similarities to that of [5] 
in the sense that it is also a hub-to-shroud inverse technique that assumes as known 

the meridional geometry and calculates the blade shape that will satisfy some flow- 
field conditions, given as input to the procedure. The input design specification used 

here is a mean swirl (radius times mean tangential velocity, rV#) distribution given 

throughout the meridional section. This design specification is somewhat unusual, but 
a suggestion in this direction can already be found in the work of Wu (see [6]) and it is 
ideally suited to the design of radial turbomachinery as discussed in [7]. In fact, the 

work of [7], which presents a three-dimensional inverse method using a mean swirl 

specification, shows that the mean swirl specification is related to the way the work is 
imparted to the fluid as it passes through the blade row. In other words, the mean swirl 
can be related to the blade loading across the blades. Another reference that discusses 

the use of a mean swirl imposition is [8], where the equations to be used in the present 
work are derived. Nevertheless [8] does not present any practical examples of 

application of the equations. 

2-DESCRIPTION OF THE DESIGN METHOD 

Throughout this work we will use a right-handed cylindrical polar coordinate 
system defined by (r, 0 , z), where r is the radius, 0 the angular coordinate and z is the 
axial distance. In addition, we will use an auxiliary coordinate a , defined by: 

a = © - f(r,z) ( 1 ) 


where f(r,z) is the angular coordinate of a point on the blade camber surface. This 
variable cc can be interpreted as a sort of helical angular coordinate aligned with the 
blade, so that when 


a = m — 
B 


( 2 ) 
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^with m an integer (m = 0, 1, 2, 3,...) and B equal to the number of blades of the 

turbomachine, we are on a blade surface (eq. (2) describes the blade shape). 

In the following we will use bold underlined characters for vectors and mean 

values will be denoted with an overbar. 

2.1 Velocity Flow Field. 

Since it was intended to apply this method to design pumps, it was assumed that 

the fluid was inviscid and incompressible and, for simplicity sake, the blade thickness 

is not considered in this procedure. 

In agreement with the through-flow approximation, the flow through the 

turbomachine will be assumed axisymmetric even in the blade region. This mean 

velocity field will be calculated using the stream function concept and the value of the 

mean vorticity. Indeed, since the vorticity field is solenoidal, it can be written as the 

cross product of two gradients of scalar functions. One of these scalar functions may 
be a according to the fact that all the vorticity is confined to the blades. In fact, if we 

suppose the far upstream velocity is uniform (an approximation quite frequent), we 

can say the flow is irrotational at inlet. Concentrating in designs that execute constant 
work along the span, it is concluded that the flow must remain everywhere 

irrotational according to Kelvin’s theorem. So, if there is any vorticity at all, it must be 
bound to the blade surfaces, justifying the statement just made. The other scalar 

function in the expression for the vorticity turns out to be the mean swirl as is shown 

in [7]. Therefore, the expression for the mean vorticity, SI is: 


SI = r Vg x V a (3) 

and now that the mean vorticity is known, the corresponding velocity field can easily 
be calculated. It is indeed known that the mean vorticity is the curl of the mean 

velocity Y . , or: 

a = vxi ( 4 ) 

Equating the 6 -component of eqs. (3) and (4) the following equation is 
obtained: 

3V r 3V Z 3f OrVfl _ 3f 3rV^ 

3z 3 r 3z 3 r 3r 3z 


relating the velocity field to the blade shape, f, and the mean swirl rV e . Besides this 
equation, the velocity field must satisfy the continuity equation. In order to achieve 
this we introduce the concept of a streamfunction defined by: 


V r 


1 3y 

r 3z 



(6a) 

(6b) 
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[so that this definition satisfies identically the continuity equation for incompressible 

flow, i. e., V . X. = 0 . The actual value of y appearing in the definitions (6a) and (6b) is 
going to be determined substituting (6a) and (6b) in eq. (5). In this way we arrive at: 

9 2 yr _ 3 2 y 

dr 2 r dr 9z 2 

For the resolution of this equation it is necessary to specify a complete set of 
boundary conditions. The boundary condition to use along the endwalls (hub and 
shroud) is the one that states that there is no flow through the solid walls. Using the 
stream function concept, this fact is expressed as: 

y = Const. (8) 

or, in other words, the hub and shroud must be streamlines of the flow. 

Far upstream we know the mean velocity vector, \L„, since it is given as input. 
Therefore we can write: 


df dr v e _ ax gr v e 

7 ?lr ^ r 


(7) 


= V 

r a s v -~ • 


IL 


(9) 


where s is distance along the far upstream boundary and n_ is the unitary vector 
perpendicular to it. This expression enables us to calculate the values of y along the 

far upstream boundary, using a simple numerical integration. At the far downstream 
boundary a similar expression applies since the velocity there is uniform because the 

flow is irrotational at inlet and the blade row is supposed to execute constant work 
along the span. In this way the complete set of boundary conditions is obtained. 

The partial differential equation (7) was solved using finite difference 
techniques. As a typical mixed-flow pump has a meridional section with complicated 

geometry bounded by curved boundaries (hub and shroud profiles) it was decided to 
use a transformation of coordinates to body-fitted curvilinear coordinates (see [9]). 
Since this coordinate system should be easy to generate and require little 
computational time, it was decided to use an algebraic transformation. For this kind of 
transformation of coordinates, mesh points are distributed along quasi-orthogonals 
and quasi- streamlines. Fig. 2 shows the grid used in the calculations to be discussed 
later on. 

Eq. (7) was discretized using second-order accurate central difference 
formulae, obtaining a nine-point difference star. The resulting finite difference 
equations were solved by a relaxation method. In the present case, a Gauss-Seidel 
relaxation scheme was used, implemented in conjunction with a multi-grid technique 
in order to accelerate the convergence rate of the solution. A good description of 
multi-grid methods can be found in [10], and, in fact, the relaxation subroutines used 
in our program are a slightly modified version of the ones presented in [10]. 

2.2 Equation for the Determination of the Blade Shape. 

After calculating the velocity field using the information presented in the 
previous subsection, it is necessary to evaluate the blade geometry. That is done by 
requiring the blade to be tangent to the velocity vector. This condition can be 
expressed as: 


YL . Vet =o 


( 10 ) 
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^where YL is the local relative velocity, =i, - II • Expanding eq. (10), the following 
expression is obtained: 



Vr dr r 2 


( 11 ) 


where f is the angular coordinate of the blade and to is the rotational speed of the blade 

row. V z and V r are the mean axial and radial velocities, respectively. 

Eq. (11) is a first-order partial differential equation with characteristic lines 
coincident with the meridional projection of the streamlines. In order to integrate this 

differential equation, some initial data must be specified along a line roughly 

perpendicular to these characteristic lines and extending from hub to shroud. This 
initial data on f will be called the stacking condition of the blade. In our method this 

stacking condition is implemented by giving, as input, the values of the blade 

coordinate f, along a quasi-orthogonal, for example, at the leading edge. 

After the stacking condition is specified, we can integrate eq. (11), since all 
the velocities that appear in it are known from the previous iteration. The integration 

of eq. (11) was done using finite difference methods. More specifically, an Euler's 

modified method (see [11]) was used. This is an implicit numerical scheme that has a 
truncation error of second order in the mesh size and is consistent and stable. 


2.3 Estimation of Blade Surface Pressures. 

One result that is important in any design method is the values of pressure on 
the suction and pressure surfaces of the blade, for judging whether the pressure 
distribution is adequate or not. In order to estimate these values starting from the 

axisymmetric solution it is necessary to calculate the blade surface velocities using the 

mean values known. To do that we begin by determining the velocity jump across the 

blades (W + - W~) which is given by (see [7]): 


w+ _ w _ _ 2n ( VrV 0 x Vq ) x Va 
B Vcc.Va 


( 12 ) 


where, for a pump, W + is the relative velocity at the pressure surface and J3L - is the 
relative velocity at the suction surface. This expression is physically plausible as it 
gives a jump in velocity which lies on the blade since it is normal to Va (a vector itself 
normal to the blade). In addition, we would expect the jump to be normal to the 


vorticity vector. 


^-[VrVfl x V a ], lying in the blade. Knowing the velocity jump and 
B 


assuming the velocity profile is linear between suction and pressure surfaces (a 
frequent approximation in hub-to-shroud methods, see [12]) it is possible to estimate 
the velocities at the suction and pressure surfaces. Indeed, their values are going to be 
equal to the mean velocity (solution of the axisymmetric problem) plus or minus one- 
half the velocity jump. Afterwards, using the fact that the flow is irrotational at inlet 
and applying Bernoulli's equation, the expression for the difference in pressure 
across the blades is obtained. The final expression is (see [7]): 


P + "P' = t p35L * VrVe 


(13) 


L 
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where p + is the pressure at the pressure surface, p~ the pressure at the suction 

surface and p is the fluid density. J5L is known from the velocity field calculation so 
that all the terms in the right-hand side of eq. (13) are known, enabling the 
estimation of the pressure loading across the blades. 

2.4 Flowchart. 

To end section 2, we would like to draw the attention of the reader to the fact 
that the vorticity depends on the blade shape, f, which in its turn is an outcome of the 
calculations. So, the calculations must be iterated until convergence is obtained. This 
is done according to the following flowchart: 

(i) input of initial data - specified values of rVg, definition of meridional section and 
body-fitted curvilinear coordinate system and all the relevant parameters; 

(ii) estimation of a first guess for f, assuming that the mean velocity is uniform along 
the quasi-orthogonals; 

(iii) using the input values of mean swirl, rVg, and the values of blade shape, f, from 
last iteration, calculate the mean velocity field (solution of eq. (7)); 

(iv) update the blade shape, f, by integration of eq. (11); 

(v) if the solution is converged output the blade shape, f, and other relevant results. 

Otherwise go back to step (iii) and initiate a new iteration. 

This flowchart was implemented as a FORTRAN computer code and applied to 
the design of the impeller of a mixed-flow pump. 

3- DISCUSSION OF APPLICATION TO A MIXED-FLOW PUMP 

In order to show the potentialities of the method, it was decided to apply it to 
the design of a mixed-flow pump. The chosen pump was based on a real machine 

which had as nominal conditions a value of 28 m for the head, H, a nominal volume 

flow, Q, of 600m^/h and a rotational speed of 1450 r.p.m.. These values give a 

nondimensional specific speed parameter equal to 0.919. The rotor has 8 blades and a 

tip diameter equal to 320 mm. This value will be used to non-dimensionalize all the 

linear dimensions and the velocities will be made nondimensional by using the 
transport blade tip velocity, wr^p (its value is 24.3 m/s). 

The meridional geometry used in the calculations is based on an existing pump 
(with minor alterations), designed by a Portuguese pump manufacturer using hand 
calculations and graphical processes. The final meridional shape used is defined in 
Fig. 1. A grid formed by 145 quasi-orthogonals and 57 quasi-streamlines was fitted to 
this meridional section, there being in the blade region a total of 61x57 points. Fig. 2 

shows every other line of the grid used. As can be seen, a region upstream and 
downstream of the blade zone was considered in the calculations. 

An important input to the present inverse method is the specification of mean 

swirl, rVg. As it is supposed that the pump accepts the flow with no swirl, the value of 

rVg along the entire leading edge is considered equal to zero. At the trailing edge the 

value was also considered constant in order to obtain a design that executes constant 

work along the span. The necessary value of rVg at the trailing edge depends on the 
work per unit mass of fluid desired for the rotor, which is a value determined by the 
desired head and an assumed value of efficiency (in our case considered equal to 0.86). 

J 


L 
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Another restriction was imposed on the derivatives of the mean swirl at the 

edge (as is requtred by .he Ku an" ze'ro presfure loading 

the pressure loading depends on 

the gradient of mean swirl, being zero where the gradient of r\ e is zero, as it is the 

CMe This ' close^ equivalence 6 between derivarives of mean swir, and pressure loading 

across the blade was one of the fac.ors used when choosing .he input mean swirl rVj, 

different turbomachine, a radial inflow turbin ^ which sufficiently 

Rprallinp hrieflv the arguments advanced in l / 1 ana wnicn <j , J 

L-jrjz mSTof *« “ 

velocity and the value of the derivative of rV„ along the meridional projection of the 
flow streamlines, or: 


+ 2k (14) 

p p ~ = T p 1 K 1 3s 

where s is distance along the meridional projection of the streamlines. In a well 

designed machine, | & | does not vary abruptly and the streamlines have a direction 
close to the quasi-streamlines. Therefore, eq. (14) implies that the pressure b 

loading is mainly influenced by the value of the derivative of rV 0 along the quasi- 

SS5VE -v “•£ srrj.'S.r.™ i 

loading with a smooth evolution, then the derivatives of r\ e should be watched with 

s«£S swSS'i !Ks 

should have o a ih ™“^ u ” r „ l i a,, ” h " at onc should bear in mind when choosing the input 

mean swirl schedule is the one connected with the i amount ^ “ 

clearly understand this argument it is important to rewrite eq. UD aog 

streamline, obtaining: 


j 
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( 15 ) 


where a and b are two arbitrary points on the same streamline, f is the angular 

coordinate of the blade, is the mean velocity in the meridional plane and m is the 
distance along the meridional projection of the streamline. If one wishes to control 
the total variation in the angular coordinate of the blade, f a - f b , it is convenient to 

avoid high values for the expression under the integral sign in eq. (15). This can be 
achieved if one specifies the value of rV# so that has small values, or in other 

words, if one specifies rVg so that the value of V# closely follows the local value of the 
transport velocity, or. This is the more important, the lower is the value of radius and 

of the meridional velocity \ m . From the point of view of highly twisted blades, the most 
critical streamsurface is the hub, not only because there the radius and meridional 

velocity take the lower values in the machine, but also because the meridional flow 
path is usually longer along the hub than anywhere else in the machine. 

The next set of four figures is presented to demonstrate that the two above 
points were taken into consideration. Indeed, Fig. 4 gives the evolution of the input 

tVq on the hub and the shroud. It is clearly seen that on the shroud the evolution of 

rV 0 is quite gradual and smooth while at the hub the opposite happens. 

Fig. 5 shows, on the same graph, the values of cor and the specified values of 

Vg at the hub. Here it is evident that, on the hub, Vq has an evolution which is almost 
parallel to cor, for most of the flow path and without much consideration in the 

direction of obtaining smooth derivatives. This was done so that at the hub 

presented small and approximately constant values along most of the blade, leading to 

a reasonable overall change in the values of blade angular coordinate, f. 

From the next figure. Fig. 6, it can be seen that the contrary happens on the 

shroud, where V# was chosen to have a smooth variation rather than following the 
local value of blade speed, cor. In fact, the main concern when specifying the mean 
swirl at the shroud was to obtain an adequate pressure distribution and not to control 

the overall variation in the blade angular coordinate, f. 

The above ideas are corroborated by the next figure, Fig. 7, where the 

derivatives of rV# along the quasi-streamlines for the hub and the shroud are 
presented. As can be clearly seen, the derivatives at the shroud present a smooth 

variation while the same does not apply at the hub, where a more abrupt change of 

the tVq derivatives can be detected. 

Using the close relationship between rV# derivatives and pressure loading one 
can conclude that the present design presents a large loading near the leading edge at 
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the shroud streamsurface, while at the hub streamsurface the blade loading is more 
evenly spread. However, one should point out that, at the hub, the derivative is zero 
for approximately the last 15% of the meridional flow path. This indicates that there 
will be a small pressure loading along the last portion of the blade at the hub, 

suggesting that perhaps one could use a shorter flow path there. This idea would entail 
an alteration of the trailing edge shape (using a trailing edge inclined to the axis, 

instead of parallel) and so, was not pursued further. 

The inverse computer code was run, using the input described above together 
with a stacking condition imposed at the trailing edge and which consists of a linear 

variation of f between the value of 0.0 (at the hub) and 0.10 (at the shroud). As a result 

the blade shape described in the next figure was obtained. Since it is difficult to 
visualize the blade three-dimensional geometry, we decided to present the blade 
geometry in Fig. 8 as a view of two consecutive blades, as would be seen by an observer 
looking in the direction of the impeller axis. The blade obtained seems typical of a 
pump impeller and no particular problems are envisaged during its manufacture, 
since it is not a highly twisted blade. 

In the next plot, Fig. 9, it is presented the estimated pressure distribution on 
hub and shroud, assuming a linear variation of the velocity from suction to pressure 
surfaces and using the procedure already discussed in subsection 2.3. The pressure 
coefficient, C p , used in this plot is defined as: 



where BL ref is a reference relative velocity, which is equal to 0.311cor tjp in the present 
case. Notice that the loading has a behaviour quite similar to the evolution of the 

derivatives of tVq along the quasi-streamlines, shown in Fig. 7, bearing out the 

comments made above, concerning the close relationship between derivatives of rVg 
and pressure loading. For example, it is seen that, at the hub, the pressure loading 

varies more abruptly than at the shroud, and along the last 15% of the flow path at 

hub the pressure loading is zero as was already expected from the values of derivatives 
of mean swirl. It should also be remarked that the distribution of pressure on the blade 
surfaces at the hub is not ideal since its variation is not smooth and presents some 
decelerations. However, the optimization of the pressure distribution would entail 
changes in the meridional section of the machine, and so was not tried in this work. 

One advantage of the present method lies in the fact that it is quite rapid, 

requiring small amounts of CPU time. In fact, the present run required lm 39s of CPU 
time in a VAX 3400 computer. Since it is computationaly so cheap, several different 
input mean swirl distributions can be scanned quickly, enabling the choice of the 
most appropriate mean swirl schedule. 

4- CONCLUSIONS 

A through-flow (hub-to-shroud) inverse method was proposed and 
implemented as a computer code. In the present implementation the flow is assumed 

incompressible, irrotational at inlet and the blade thickness was not considered 

during the calculations. The necessary equations are presented and, as an example of a 
possible application, the inverse method was used to redesign the rotor of a mixed-flow 
pump. 

The present method uses as an input specification the value of mean swirl, 

rVfl . This input specification was chosen using a reasoning similar to that used in [7], 
This is remarkable, since the turbomachine designed in [7] was a radial-inflow 
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turbine, and indicates that the strategy developed there is quite general, being able to 
cope with radically different machines. Briefly, this strategy consists in using the 
input mean swirl to control the pressure loading on the blades and the overall 
variation in the angular coordinate of the blade, f. 

This work also shows that further research should be done in order to obtain 
reasonable pressure distributions on the blade surfaces. In fact, the design presented 
here has a pressure distribution on the hub which is not the ideal one. The 
optimization of the pressure distribution will probably involve some changes in the 
meridional section (hub and shroud contours). 

One advantage of the present technique is the fact that it requires small 
amounts of CPU time. So it is a convenient tool to scan quickly and inexpensively 
several different input mean swirl distributions, in order to find the most appropriate 
one. 
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